N 7 2 “ 23 195 


NASA TECHNICAL 

MEMORANDUM 


NASA TM X-68041 


iH 

s 

oo 

vO 


X 






ILE 


PILOT/VEHICLE CONTROL OPTIMIZATION USING AVERAGED 
OPERATIONAL MODE AND SUBSYSTEM RELATIVE 
PERFORMANCE INDEX SENSITIVITIES 


by Gary G. Leininger 
University of Toledo 
Toledo, Ohio 

and 

Bruce Lehtinen, John P. Riehl 
Lewis Research Center 
Cleveland, Ohio 

TECHNICAL PAPER proposed for presentation at 
Eighth Annual National Aeronautics and Space 
Administration - University Conference on Manual Control 
Ann Arbor, Michigan, May 17-19, 1972 


E-6879 


PILOT/VEHICLE CONTROL OPTIMIZATION USING AVERAGED 
OPERATIONAL MODE AND SUBSYSTEM RELATIVE 

PERFORMANCE INDEX SENSITIVITIES 

\ 

by Gary G. Leininger, * Bruce Lehtinen, ^ and John P. Riehl^ 

ABSTRACT 

A method is presented for designing optimal feedback controllers 
for systems having subsystem sensitivity constraints. Such constraints 
reflect the presence of subsystem performance indices which are in con- 
flict with the performance index of the overall system. The key to the 
approach is the use of relative performance index sensitivity (a measure 
of the deviation of a performance index from its optimum value). The 
weighted sum of subsystem and/or operational mode relative perform- 
ance index sensitivies is defined as an overall performance index. A 
method is developed to handle linear systems with quadratic perform- 
ance indices and either full or partial state feedback. The usefulness 
of this method is demonstrated by applying it to the design of a stability 
augmentation system (SAS) for a VTOL aircraft. A desirable VTOL SAS 
design is one that produces good VTOL transient response both with and 
without active pilot control. The system designed using the method in- 
troduced in this paper is shown to effect a satisfactory compromise solu- 
tion to this problem. 
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1.0 INTRODUCTION 

The study of complex systems often involves an investigation into 
the interconnection of many subsystems and the influence each subsys- 
tem has in achieving a prespecified design objective. The optimization 
of the composite system with respect to a set of adjustable parameters 
relies upon a knowledge and understanding of the interconnecting struc- 
ture. Individual subsystem optimization without concern for its cause 
and effect relation to the composite system may yield an overall sys- 
tem response which deviates substantially from the design specifica- 
tions. Conversely, a composite system design satisfying the required 
design objectives might dictate a need for increasingly sophisticated and 
expensive subsystems. It is, therefore, in the interests of the system 
designer to have the practical and analytical flexibility to properly 
aline the priorities of competitive subsystems with composite system 
objectives. A system typifying this design analysis is a pilot/vehicle 
system. 

In the process of designing and evaluating the suitability of a 
pilot/vehicle system, it is necessary to solicit the pilot's comments 
and opinion of the handling qualities as one facet of the design. This 
subjective opinion forms an integral part of the ultimate evaluation of 
the vehicle and is therefore considered seriously and continuously 
.throughout the design^ ‘Optimal performance- of the-airc-raft--including 
the pilot may be in direct conflict with the optimal performance and ef- 
ficiency of the aircraft when the aircraft is treated as a separate en- 
tity. In this respect many studies have been conducted with the objec- 
tive of mathematically modeling pilot-control characteristics (refs. 1 
to 3) and from a practical engineering viewpoint, the development of 
the quasi-linear model for human pilot dynamics have been one of the 
beneficial results of these studies (ref. 1). These pilot models can 
then be used in conjunction with airframe dynamic models in the design 
of aircraft control systems. 
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Composite system design can best be achieved when the design 
criterion includes the evaluations, requirements, and limitations of 
each individual subsystem. To accomplish this, a generalized theory 
and design technique is presented. The theory evolves from the con- 
cepts and conditions imposed by optimal control theory supplemented by 
subsystem sensitivity characteristics. 

The objective of optimal control theory is to determine the control 
signals that will cause a process to satisfy the physical constraints 
and at the same time minimize some performance criterion. In the 
case of feedback control, the parameters to be optimized are the feed- 
back gains. Once the optimization has been completed, it is natural 
to inquire into the relative effect of the system response and/or the 
performance measure to a deviation of the feedback gains from their 
optimal values. This area of concern is often referred to as "sensi- 
tivity. " 

The design and evaluation of dynamic control systems through the 
utilization of sensitivity functions has been the subject of intensive 
research during the past decade. Many different definitions of sensi- 
tivity have evolved and system stability, controllability, and other 
system characteristics have been directly related to these sensitivity 
functions. 

Many analysis, synthesis and optimization techniques used in con- 
trol theory utilize the sensitivity functions of the state of the system 
with respect to the system parameters (refs. 4 and 5). These param- 
eter sensitivity functions are often generated by sensitivity models of 
the system. However, the use of presently available techniques for 
generating sensitivity functions for linear system containing many pa- 
rameters results in the simulation of high order dynamic systems. 
Similarly, control sensitivity (ref. 6), trajectory sensitivity (refs. 

4 and 7), eigenvalue sensitivity (ref. 8) and output sensitivity with re- 
spect to pole location (ref. 9) all require a simultaneous solution of 
high order dynamic sensitivity expressions. These sensitivity meth- 
ods, however, are only intended to account for very small perturbations 
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from some nominal (or optimal) position. For this reason, most of the 
above sensitivity functions are not compatible with on-line design, but 
are primarily used to evaluate the final design. Papers by Cadzow 
(ref. 10), Dougherty, et al. (ref. 11) and others have applied perform- 
ance index sensitivity methods to the problem of determining feedback 
control laws when system parameters are subject to small variations. 
For larger variations, authors such as Whitbeck (ref. 12), Zadicario 
and Sivan (ref. 13), and Tuel (ref. 14) have discussed methods for de- 
signing controllers which minimize the expected value of the cost func- 
tional. All of these techniques are addressed to the plant parameter 
variation problem as well as a single scalar performance measure. 

The emphasis in this paper is directed toward developing a method 
for designing a practical feedback controller for multivariable linear 
systems which may be stabilized by output feedback over the entire 
range of feedback parameters. The distinguishing feature of the pro- 
posed technique is the generation of a constant feedback control law 
subject to the minimization of the performance index sensitivity func- 
tions of the composite system and the individual subsystems. This is 
accomplished by defining a performance index consisting of the sum of 
relative sensitivity terms of each subsystem multiplied by scalar 
weighting factors. 

Relative sensitivity is a measure of the deviation between the 
actual value of the performance index and that which would be obtained 
if the control were optimal, i. e. , 


s R(k) a m - j(K° i 1 (1 .D 

J(K°) 

where K° is the optimal set of feedback parameters 

J(K°) = min J(K) (1.2) 

K° 
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Note that the relative sensitivity is always positive, and thus system 
performance is always compared with an attainable value. 

The performance index of concern here is chosen to be of the form 

N. 

J(K)= £ XjSf(K) (1.3) 

i=l 

t 

N 

Yj \ - 1 (1. 4) 

i=l 

where N is the number of subsystems, S^(K) is the relative sensitivity 
of the i un system, and X. is a weighting factor (or probability factor) 
associated with the i subsystem. The performance index of equa- 
tion (1. 3) reflects the interest and concerns of the individual subsystems 
in the overall decision process. Clearly, small relative sensitivity 
assures a design close to the optimum and, hence, a smaller influence 
in the final optimization procedure. Furthermore, the performance 
measure J^K) associated with subsystem i (i = 1, 2, . . . , N) need 
not be of the same form, i. e. , quadratic, absolute value, uniform, etc. 
This, therefore, greatly enhances the design capabilities for large scale 
systems with subsystem design limitations. 

The merit of the performance index of equations (1. 3) and (1. 4) in 
the design of practical engineering systems is considered in the next 
section. Here we restrict the discussion to linear systems with quad- 
ratic performance criteria. This restriction enables the designer to 
utilize the well-developed theory of the optimal linear regulator in es- 
tablishing the sensitivity terms needed in equation (1. 3). The section is 
then concluded with an example demonstrating the effectiveness of the 
technique as a useful design tool. 
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2. 0 SENSITIVITY DESIGN FOR OPTIMAL 
LINEAR REGULATORS 

This investigation is concerned with the design of a feedback con- 
trol law for a time invariant linear system subject to the minimization 
of a prespecified scalar performance index. The performance index 
is chosen in such a manner as to include sensitivity terms associated 
with subsystems comprising the composite system. Initial considera- 
tion is given to the optimization of an individual subsystem with respect 
to a selected array of feedback parameters. The method is then ex- 
tended to include several subsystems in the overall optimization proce- 
dure using sensitivity concepts derived for the individual subsystem. 

Subsystem Optimization 

Consider a subsystem whose dynamic performance is character- 
ized by a set of n first order linear time invariant differential equa- 
tions. 

x(t) = Ax(t) + Bu(t) (2. 1) 

•where x(t) is the n dimensional state -vector, u(t)=isa vector consist- 
ing of m control inputs, and A and B are nxn and nxm constant 
matrices describing the system dynamics. The feedback control law 

u(t) = -Lx(t) (2.2) 

will be optimal if the feedback gain matrix L is chosen so as to mini- 
mize a performance index which is quadratic in the state and control 
variables 
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(2. 3) 

where weighting matrix Q is positive semidefinite and weighting ma- 
trix R is positive definite. The minimization of equation (2. 3) will 
yield a set of constant feedback gains of the form 

L* - R -1 B T P* (2.4) 

where P* is a positive definite matrix which is the solution of the 
steady- state Riccati equation 

0 = A T P* + P*A + Q - P*BR -1 B T P* (2. 5) 

■s 

Combining equations (2. 4) and (2. 5), one obtains 

0 “ (A - BL*) T P* + P*(A - BL*) + Q + L* T RL* (2. 6) 

which is the well-known Lyapunov equation. The resulting value of 
the performance index when the feedback gain matrix equation (2. 4) is 
substituted into equation “(2. 2)' is * ‘ 

J *=^ X 0 P * X 0 (2 - 7) 

The evaluation of P* from equation (2. 5) requires the solution of 
n(n + l)/2 nonlinear simultaneous equations. Alternatively, one could 
use equations (2. 6) and (2. 7) in conjunction with a gradient minimiza- 
tion algorithm. For a nonoptimal set of feedback gains L, equa- 
tion (2. 6) becomes 

0 = (A - BL) T P + P(A - BL) + Q + L T RL 


J(x, u) — — 
2 


J 0 


{x T (t)Qx(t) + u T (t)Ru(t)}dt 


( 2 . 8 ) 
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Here A, B, Q, R, and L are known and the nxn symmetric matrix P 
can be easily obtained using any of the well-known Lyapunov solving 
algorithms (refs. 16 and 19). The performance index is evaluated using 

J=IxT P x o (2.9) 

and can be minimized by adjusting the elements of the gain matrix using 
a gradient minimization algorithm yielding L* and P*. 

In most practical situations, the initial state of the system is un- 
known and must be treated as a random vector. Taking the expected 
value of equation (2. 9) yields 


n 

E { J } = 3 = £ y u p u (2.10) 

i=l 

where p.. are the diagonal elements of the P matrix and y .. is the 

XL in X X 

covariance of the i L component of the initial state vector with the 
additional assumption that 

Vy =E<x 1 (0)x j (0)} =0 (2.11) 

Note that the minimization of equation (2. 10) subject to equation (2. 8) 
will yield an 5 ’averaged” set of feedback gains independent of the statis- 
tics of the initial state random vector. 

The above analysis assumes that the feedback control law u(t) is a 
linear combination of all the elements of the state vector. In the event 
that only a select number of the state variables will comprise the feed- 
back control law, equation (2. 8) must be modified. Let the p dimen- 
sional vector 


y(t) = Cx(t) 


( 2 . 12 ) 
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represent the state variables to be fed back. For a control law of the 
form 


u(t) = -Ky(t) = -KCx(t) (2. 13) 

equation (2. 8) becomes 


0 = (A - BKC) T P + P(A - BKC) + Q + C T K T RKC (2. 14) 

The minimization of equation (2. 10) subject to equation (2. 14) proceeds 
as above with the additional restriction that the closed loop system 
(A - BKC) be stable. This latter restriction is, of course, in effect in 
the full state feedback system; however, it is well known that the op- 
timal linear state regulator is always stable independent of the open 
loop dynamics. Clearly, this is not true, in general, for the partial 
state feedback system and thus one must be cognizant of the location of 
the closed loop poles, since any solution to equation (2. 4) yielding an 
unstable closed loop system is meaningless. 

Relations similar to those of equations (2. 8) and (2. 14) have been 
obtained by Kleinman (ref. 17) and Levine (ref. 18), respectively. 
However, their results are predicated upon the existence of the first 
partial derivative of equation .(2,. 10) with respect.to the .unknown feed- 
back gain matrix. Setting this derivative to zero provides the relations 
upon which their derivation and subsequent results ultimately rely. 
Consequently, if the feedback gains are constrained in any manner, 
then 3J/3L + 0 or 3J/3K + 0 at the optimum and the results of 
(ref. 17) and (ref. 18) no longer apply. However, the constraint 
boundaries can be incorporated into the gradient algorithm described 
above and thus minimization of equation (2. 10) with respect to the 
constrained gains can be achieved. The stabilization of (A - BL) for 
full state feedback and (A - BKC) for partial state feedback is an addi- 
tional restriction placed upon the feedback gains. 
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Up to this point a technique has been presented for the optimiza- 
tion of a single subsystem with respect to a set of feedback gains com- 
prising a feedback control law. The method was shown applicable to 
systems with full state feedback, partial state feedback and both full 
and partial state feedback with gain constraints. The method will now 
be extended to include several subsystems with conflicting objectives. 

Consider a composite system ^“q which contains a definable sub- 
system of interest cTg- Two quadratic performance indices, Jq and 
J 2 , are defined for c^q and cT 2 , respectively. Each index has been 
suitably averaged over the initial conditions (as was done in eq. (2. 10)). 
For ease of discussion, assume the system is structured as follows 
(see figs. 1 and 2). 


System gTq 


A 1 A 
A 11 1 A 12 


0 


1 A 
1 A 22 


X 1 

I---1 + 


_0 

B 


'n 


1 l 2 


0 

B, 


' u - 


vu« 


y_i_ 

\ y 2, 


Ci 0 


0 


|C 2 


Xr 


u< 


u. 


-K 1 1 0 
1 l 

~0~ '~Kr 


72 . 


.(2.-15). 


min J Q (x 2 , Up u 2 ) J j* 


k p k 2 


, / 

\ 

\: 
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Subsystem c/T^g (System c/°g with 0^ = 0) 


x 2 “ A 22 x 2 + B 2 u 2 


y 2 = C 2 X 2 


u 


2 " “ K 2 y 2 


A A 

min J g(xg, u 2^ = J 2 
K n 


} ( 2 . 16 ) 


Note that both Jg and Jg are functions of subsystem states (X 2 ) only. 
0 0 

Let K^ and Kg be the values of and Kg that minimize Jg, and 

2 - , 

let Kg be the value of Kg that minimizes Jg. In most physical 

situations Kg ± Kg. This fact reflects a degradation in the perform- 
ance of subsystem cTg when the gains Kg that optimize system ci?g 
are used. This degradation can be measured using relative sensitivity, 
i. e. , 



( 2 . 17 ) 


where Jg is the value of the performance index for subsystem c/°g 
when the feedback gains Kg are employed. A similar expression for 

a, \ 

the relative sensitivity of Jg is 


S^(K) 




( 2 . 18 ) 
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where Jq is the value of Jq at Kp K^. Clearly, = 0 if 

K 2 = K 2’ an£ * Sq ~ 0 if Kj = . and = K®. Consequently, a 

composite performance index can be defined which incorporates equa- 
tions (2. 17) and (2. 18). 

J 3 = XS^(K) + (1 - X)S**(K) (2. 19) 

where X is a weighting (or probability) factor 0 < X < 1. For X = 0, 
the minimization of equation ( 2 . 18) will yield the gains and K° 
while for X = 1 the gains will result and will have no effect 
in the minimization. Thus, for X in the range of zero to one the 
minimization of equation (2. 18) will result in a tradeoff between the 
design objectives of subsystem and system c?q. The appropriate 
value of X will depend entirely upon the physical, as well as design, ' 
requirements of both the individual subsystem and the composite sys- 
tem. 

To generalize the above technique, consider a composite system 
<£q, and N definable subsystems cSX, i — 1, 2, . . . , N, which may 
be coupled either through state or through control. For each subsys- 
tem, define a relative sensitivity of the form 

R J i " J t 

■ ' SfXK) W V 1 - ‘ * (2. 20) 

J i 

where J* is the value of the performance index when (sub)system i is 
optimized independent of the rest of the system. The relative perform- 
ance index can then be formulated as a linear combination of the sub- 
system sensitivity functions 

j - £ \ s f< K > 

i«0 


( 2 . 21 ) 
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N 

Y) x i = 1 (2. 22) 

i=0 

th 

where is the weighting factor associated with the i in (sub)system. 
The minimization of equation (2, 21) will yield a set of feedback gains 
for the composite system which will provide satisfactory overall per- 
formance, while maintaining subsystem response within the design (and 
economic) specifications,, To demonstrate the effectiveness of the pro- 
posed technique as a suitable design tool, as well as to illustrate the 
approach, an example is now presented. 

3.0 EXAMPLE 

Consider the problem of designing the stability augmentation sys- 
tem (SAS) for a turbojet/lift fan powered VTOL aircraft. Figure 3 
is the block diagram for a linearized pitch axis model of a typical 
VTOL being controlled by a pilot in the hover mode. The states con- 
sidered in the VTOL model are pitch angle, pitch rate, and accelera- 
tion produced by the moment generated by the lift fans. To be designed 
are SAS gains K^, and such that pitch angle 9 is kept close 
to zero using reasonable amounts of control, u. The complete system 
is similar in form to the one depicted in figure 2, where the VTOL 
pitch dynamics comprise the primary subsystem of interest. In this 
example, reflects the fact that only one pilot state (6) is measur- 
able; and for simplicity, the gain associated with this state was 
fixed. 

Pilot dynamics are described by the third order model shown in 
the diagram (ref. 20). Pilot parameters are pilot gain Kp, lead time 
constant r^, muscle lag Tjyp sensor lag Tg, and pilot dead time 
(as a Pade approximation). Aircraft handling qualitites studies have 
shown that if the parameters the pilot adaptively adjusts (Kp and t l ) 
are not too large, the pilot will give the aircraft a high rating. For 
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this example, Kp was fixed at 13. 5 cm/rad and r L was fixed at 
0. 5 sec. These are relatively low values; hence the resulting design 
should get a good pilot rating. Fixed pilot parameters were assumed 
to be Tg = 0. 062 sec, r M - 0. 36 sec, and r D = 0. 35 sec. As indi- 
cated previously, stick sensitivity K* was not optimized, but was 

2 ^ 

chosen to be 0„ 6 rad/sec /in. , based on typical pilot preferences. 

The VTOL dynamic model includes parameter which repre- 
sents the lag between stick deflection and pitching moment produced 
by the engine/lift fan combination, j ^ was assumed to be 0. 3 second 
in this study, since experience has shown that pilot may have diffi- 
culty in controlling the system is the actuation lag is greater than this 
amount. Conversely, it is desirable, as far as lift fan/engine design 
is concerned, to have r- ^ as large as possible. 

The VTOL can be operated in either of three modes: (1) pilot-in- 
the-loop (PIL) where both pilot and SAS contribute to stabilization 
(Kp * 0), and (2) pilot-out-of-the-loop (POL) where Kp = 0 and all 
stabilization derives from the SAS, and (3) SAS-failed mode, where 
the pilot provides all stabilization. Because of these three modes of 
possible operation, a conflict arises in designing the SAS and stick 
sensitivity K g . It was decided that K g would not be optimized in 
this study, hence, the SAS-failed mode has been ignored. For the PIL 
and POL modes, SAS optimized for the POL mode may produce a sys- 
tem too insensitive to the pilot's control during PIL mode operation. 

On the other hand, a SAS which is designed to be optimal when the pilot 
is in the loop may not sufficiently stabilize the aircraft in the fixed- 
stick (POL) mode. For this example, the problem of conflicting per- 
formance objectives was solved by using a composite performance index 



(3.1) 


A /S 

Here, J 2 is an index of performance for POL operation, Jq corre- 
sponds. to PIL operation, and \ is a weighting factor. Lambda could 
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be, for example, the probability of the aircraft being flown in the POL 
mode. The starred J values are those that are obtained when optimiz- 
ing for PIL or POL mode operation separately. 

The form of the performance index chosen to be minimized in each 
mode of operation is 


J a El J} = E 


In 


9 2 + k 6 6 2 + kjjU^g Idt 


(3. 2) 


where k c 


and k u are scalar weighting constants and J is to be 


averaged over the initial states. The three terms in this performance 
index were selected in accordance with the following considerations: 

(1) d should be driven to zero as rapidly as possible, (2) required pilot 
stick deflection should not be excessive, (3) a control moment command, 
UgAg, generated by the SAS should not cause the lift fans to exceed 
their rated thrust. For POL operation, the aircraft performance index 
can be written in the form of equation (2. 10) as 


J 2 - Yj X ii P ii 

i=l 


(3. 3) 


where the P matrix for the aircraft subsystem 'is obtained "By solving 
a third order equation of the form of equation (2. 14). Similarly, for 
the pilot-in-the-loop, 

J o= EniPu (3 - 4) 

i=l 

Here P is the solution of equation (2. 14) for the complete sixth order 
system with the stipulation that y.. “ 0 for the three pilot states 
(initial pilot states are assumed to have zero mean and variance). For 
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both modes, initial states are assumed to be uncorrelated (y. . = 0 , 

i*j)- 

Solutions were obtained (optimal SAS gains) for various values 
of X. Powell's method (ref. 21) of function minimization was used 
along with the Lyapunov equation solution technique of reference 16 for 
evaluation J. Covariances y 722 ’ an( * V 33 were a ll assumed to 
be 1,0, and weighting factors kg and k u were chosen as 0. 0015 
and 0„ 15, respectively. Figure 4 graphically presents the results of 
the optimization of the composite performance index, J^, as a func- 
tion of SAS pitch gain for a selection of X values. Each curve 
is a section through the performance surface with and Kg held 
constant at the optimum values obtained for that particular X value. 
Optimal Kg s (which occur at the minima) range from 2. 19 to 3. 69; 
however, their magnitudes are quite similar for the extreme cases 
(X = 0 and 1). The design trade-off is evidenced by the fact that 
increases as X moves away from 0 or 1, up to a maximum of about 
0. 04 for the X = 0. 6 curve. 


To further demonstrate the influence of X on control system 
behavior, typical transient responses for fixed initial conditions were 
computed. Figures 5(a) and (b) show, respectively, PIL and POL 
responses of pitch angle 9 , moment command due to the pilot u , 

r 

moment command due to the SAS Ug^g, and the resulting VTOL pitch 
acceleration 9 . Transients are displayed for four X values (0, 0. 2, 

0. 8 , and 1. 0) for zero pilot initial conditions and VTOL initial condi- 
tions of 0 . 1 , 0 . 1 , and 0 . 1 . 

Comparing the pitch angle transients, it can be seen that for PIL 
operation (fig. 5(a)), the best transient occurs for the system optimized 
for X = 0. Conversely, the best POL transient (fastest response) 
occurs in figure 5(a) for the system optimized for X = 1. As an ex- 
ample, consider the case when the VTOL is in the POL mode 80 percent 
of the time, i. e. , X = 0. 8 . The 9 curve for X = 0. 8 in figure 5(a) 
shows performance is somewhat degraded (higher overshoot, poorer 
damping) over the X = 0 case, , but is not nearly so poor as the highly 
underdamped X = 1 case. In figure 5(b) it can be seen that the X = 0. 8 . 
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design is very nearly as fast responding as the X - 1 design, and much 
faster than the system designed for PIL operation (X =0). 


Similar comparisons can be made for u , u s^S’ an< * ® transients. 
For instance, pilot control excursion is high in figure 5(a) for the X = 1 


case, (PIL operation with POL feedback gains). However, using these 


same gains in POL operation gives the best transient performance. 


Note that Up is zero in figure 5(b), since the pilot is not exercising 
control. Pitch acceleration histories are included to demonstrate that 


all of the optimal controllers give rise to VTOL accelerations which 
are '’reasonable’ 8 in magnitude. The rather anomalous behavior of Up 
in PIL operation (fig. 5(a)) at t - 0 is due to the fact that in the analysis 
and transient calculation, pilot dead time has been approximated by a 
first order Pade approximation. What appears to be the pilot initially 
attempting to increase the error in 0 is actually due to the inaccuracy 
in modeling his dead time with a Pade. 

In the preceeding example, no constraints were imposed on the SAS 
gains. One obvious constraint that could be considered is one on Kg. 

As becomes large, the lift fan/engine eigenvalue increases, such 
that eventually saturation will certainly occur. Thus, for a reasonable 
solution, must be bounded. Another problem, mentioned in Sec- 
tion 2, is system stability. The system in this example was open loop 
stable such that even though not all states were fed back, a set of (op- 
timal) gains were found which produced a stable system. This will not 
be the case, in general, so that periodic stability checks must be made 
during the optimization to insure each set of (sub-optimal) gains corre- 
sponds to a stable system. 


4.0 CONCLUSIONS 

An approach has been formulated to the problem of designing a con- 
trol for a system with conflicting subsystem performance indices. Use 
was made of relative sensitivity by introducing it into the system’s per- 
formance index. A method was developed for handling linear systems 



18 


with quadratic subsystem performance indices, for either full or 
partial state feedback. The approach was demonstrated by using it 
to design the pitch axis SAS for a piloted VTOL, where the main sub- 
system of interest was the VTOL aircraft. A design was obtained, 
consisting of a fixed set of SAS gains, which gave acceptable perform- 
ance both with and without pilot control. The methods developed could 
be extended to include nonlinear plants, state variable constraints, 
and nonquadratic performance indices. They could also be applied to 
designing the complete three axis SAS for a VTOL, capable of operat- 
ing throughout the hover, transition and cruise modes. 
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